ALK-rearranged renal cell carcinoma with TPM3::ALK gene fusion and review of the literature.

ALK-rearranged renal cell carcinoma (ALK-RCC) is a very rare novel molecularly defined entity in the recently published fifth edition of the World Health Organization classification of tumours. We describe a case of ALK-RCC in a 76-year-old female. The tumour was composed of discohesive rhabdoid cells and pleomorphic, multinucleated cells (equivalent to ISUP/WHO grade 4). The tumour showed expression with PAX8, Keratin 7 and alpha methylacyl CoA racemase. ALK (D5F3 clone) was strongly and diffusely positive. ALK-FISH showed significant split signals of ALK, confirming the diagnosis. RNA sequencing showed TPM3::ALK rearrangement. Including the current case, there are 14 reported ALK-RCC cases with the same TPM3 fusion partner gene. Review of these published cases highlights their morphological heterogeneity and stresses the importance of running ALK immunohistochemistry on difficult cases to classify renal tumours. This is important while identification of ALK-RCC has clinical significance due to the availability of targeted therapy with ALK inhibitors.

A MXI1-NUTM1 fusion protein with MYC-like activity suggests a novel oncogenic mechanism in a subset of NUTM1-rearranged tumors.

Most NUTM1-rearranged neoplasms (NRNs) have fusions between NUTM1 and BRD (bromodomain-containing) family members and are termed NUT carcinomas (NCs) because they show some squamous differentiation. However, some NRNs are associated with fusions between NUTM1 and members of the MAD (MAX dimerization) gene family of MYC antagonists. Here we describe a small round cell malignancy from the gastro-esophageal junction with a previously unreported fusion between NUTM1 and the MAD family member MXI1. In contrast to NCs, the MXI1-NUTM1 tumor did not show squamous differentiation and did not express MYC, TP63 or SOX2, genes known to be targets of BRD-NUTM1 proteins and critical for NC oncogenesis. Transcriptome analysis showed paradoxical enrichment of MYC target genes in the MXI1-NUTM1 tumor despite the lack of MYC expression. When expressed in vitro MXI1-NUTM1 partially phenocopied MYC, enhancing cell proliferation and cooperating with oncogenic HRAS to produce anchorage-independent cell growth. These data provide evidence that MAD family members, which are normally repressors of MYC activity, can be converted into MYC-like mimics by fusion to NUTM1. The pathological features and novel oncogenic mechanism of the MXI1-NUTM1 tumor show that identification of NUTM1 fusion partners can be important for accurate diagnostic classification of some NRN subtypes, and potentially may guide therapeutic options.

RAF1 rearrangements are common in pancreatic acinar cell carcinomas.

There is now evidence that gene fusions activating the MAPK pathway are relatively common in pancreatic acinar cell carcinoma with potentially actionable BRAF or RET fusions being found in ~30%. We sought to investigate the incidence of RAF1 fusions in pancreatic malignancies with acinar cell differentiation. FISH testing for RAF1 was undertaken on 30 tumors comprising 25 'pure' acinar cell carcinomas, 2 mixed pancreatic acinar-neuroendocrine carcinomas, 1 mixed acinar cell-low grade neuroendocrine tumor and 2 pancreatoblastomas. RAF1 rearrangements were identified in 5 cases and confirmed by DNA and RNA sequencing to represent oncogenic fusions (GATM-RAF1, GOLGA4-RAF1, PDZRN3-RAF1, HERPUD1-RAF1 and TRIM33-RAF1) and to be mutually exclusive with BRAF and RET fusions, as well as KRAS mutations. Large genome-wide copy number changes were common and included 1q gain and/or 1p loss in all five RAF1 FISH-positive acinar cell carcinomas. RAF1 expression by immunohistochemistry was found in 3 of 5 (60%) of fusion-positive cases and no FISH-negative cases. Phospho-ERK1/2 expression was found in 4 of 5 RAF1-fusion-positive cases. Expression of both RAF1 and phospho-ERK1/2 was heterogeneous and often only detected at the tumor-stroma interface, thus limiting their clinical utility. We conclude that RAF1 gene rearrangements are relatively common in pancreatic acinar cell carcinomas (14.3% to 18.5% of cases) and can be effectively identified by FISH with follow up molecular testing. The combined results of several studies now indicate that BRAF, RET or RAF1 fusions occur in between one third and one-half of these tumors but are extremely rare in other pancreatic malignancies. As these fusions are potentially actionable with currently available therapies, a strong argument can be made to perform FISH or molecular testing on all pancreatic acinar cell carcinomas.

Emerging entities in NUTM1-rearranged neoplasms.

Structural alterations of NUTM1 were originally thought to be restricted to poorly differentiated carcinomas with variable squamous differentiation originating in the midline organs of children and adolescents. Termed NUT carcinomas (NCs), they were defined by a t(15;19) chromosomal rearrangement that was found to result in a BRD4-NUTM1 gene fusion. However, the use of DNA and RNA-based next-generation sequencing has recently revealed a multitude of new NUTM1 fusion partners in a diverse array of neoplasms including sarcoma-like tumors, poromas, and acute lymphoblastic leukemias (ALLs) that we propose to call NUTM1-rearranged neoplasms (NRNs). Intriguingly, the nosology of NRNs often correlates with the functional classification of the fusion partner, suggesting different oncogenic mechanisms within each NRN division. Indeed, whereas NCs are characterized by their aggressiveness and intransigence to standard therapeutic measures, the more positive clinical outcomes seen in some sarcoma and ALL NRNs may reflect these mechanistic differences. Here we provide a broad overview of the molecular, nosological, and clinical features in these newly discovered neoplastic entities. We describe how aberrant expression of NUTM1 due to fusion with an N-terminal DNA/chromatin-binding protein can generate a potentially powerful chromatin modifier that can give rise to oncogenic transformation in numerous cellular contexts. We also conclude that classification, clinical behavior, and therapeutic options may be best defined by the NUTM1 fusion partner rather than by tumor morphology or immunohistochemical profile.

Improved next-generation sequencing pre-capture library yields and sequencing parameters using on-bead PCR.

Tumor DNA sequencing results can have important clinical implications. However, its use is often limited by low DNA input, owing to small tumor biopsy size. To help overcome this limitation we have developed a simple improvement to a commonly used next-generation sequencing (NGS) capture-based library preparation method using formalin-fixed paraffin-embedded-derived tumor DNA. By using on-bead PCR for pre-capture library generation we show that library yields are dramatically increased, resulting in decreased sample failure rates. Improved yields allowed for a reduction in PCR cycles, which translated into improved sequencing parameters without affecting variant calling. This methodology should be applicable to any NGS system in which input DNA is a limiting factor.

Development and interlaboratory evaluation of a NIST Reference Material RM 8366 for EGFR and MET gene copy number measurements.

Background The National Institute of Standards and Technology (NIST) Reference Material RM 8366 was developed to improve the quality of gene copy measurements of EGFR (epidermal growth factor receptor) and MET (proto-oncogene, receptor tyrosine kinase), important targets for cancer diagnostics and treatment. The reference material is composed of genomic DNA prepared from six human cancer cell lines with different levels of amplification of the target genes. Methods The reference values for the ratios of the EGFR and MET gene copy numbers to the copy numbers of reference genes were measured using digital PCR. The digital PCR measurements were confirmed by two additional laboratories. The samples were also characterized using Next Generation Sequencing (NGS) methods including whole genome sequencing (WGS) at three levels of coverage (approximately 1 ×, 5 × and greater than 30 ×), whole exome sequencing (WES), and two different pan-cancer gene panels. The WES data were analyzed using three different bioinformatic algorithms. Results The certified values (digital PCR) for EGFR and MET were in good agreement (within 20%) with the values obtained from the different NGS methods and algorithms for five of the six components; one component had lower NGS values. Conclusions This study shows that NIST RM 8366 is a valuable reference material to evaluate the performance of assays that assess EGFR and MET gene copy number measurements.

Profound MEK inhibitor response in a cutaneous melanoma harboring a GOLGA4-RAF1 fusion.

BRAF and CRAF are critical components of the MAPK signaling pathway which is activated in many cancer types. In approximately 1% of melanomas, BRAF or CRAF are activated through structural arrangements. We describe here a metastatic melanoma with a GOLGA4-RAF1 fusion and pathogenic variants in CTNNB1 and CDKN2A. Anti-CTLA4/anti-PD1 combination immunotherapy failed to control tumor progression. In the absence of other actionable variants the patient was administered MEK inhibitor therapy on the basis of its potential action against RAF1 fusions. This resulted in a profound and clinically significant response. We demonstrated that GOLGA4-RAF1 expression was associated with ERK activation, elevated expression of the RAS/RAF downstream co-effector ETV5, and a high Ki67 index. These findings provide a rationale for the dramatic response to targeted therapy. This study shows that thorough molecular characterization of treatment-resistant cancers can identify therapeutic targets and personalize management, leading to improved patient outcomes.

Breast cancer metastasis to gynaecological organs: a clinico-pathological and molecular profiling study.

Breast cancer metastasis to gynaecological organs is an understudied pattern of tumour spread. We explored clinico-pathological and molecular features of these metastases to better understand whether this pattern of dissemination is organotropic or a consequence of wider metastatic dissemination. Primary and metastatic tumours from 54 breast cancer patients with gynaecological metastases were analysed using immunohistochemistry, DNA copy-number profiling, and targeted sequencing of 386 cancer-related genes. The median age of primary tumour diagnosis amongst patients with gynaecological metastases was significantly younger compared to a general breast cancer population (46.5 versus 60 years; p < 0.0001). Median age at metastatic diagnosis was 54.4, time to progression was 4.8 years (range 0-20 years), and survival following a diagnosis of metastasis was 1.95 years (range 0-18 years). Patients had an average of five involved sites (most frequently ovary, fallopian tube, omentum/peritoneum), with fewer instances of spread to the lungs, liver, or brain. Invasive lobular histology and luminal A-like phenotype were over-represented in this group (42.8 and 87.5%, respectively) and most patients had involved axillary lymph nodes (p < 0.001). Primary tumours frequently co-expressed oestrogen receptor cofactors (GATA3, FOXA1) and harboured amplifications at 8p12, 8q24, and 11q13. In terms of phenotype conversion, oestrogen receptor status was generally maintained in metastases, FOXA1 increased, and expression of progesterone receptor, androgen receptor, and GATA3 decreased. ESR1 and novel AR mutations were identified. Metastasis to gynaecological organs is a complication frequently affecting young women with invasive lobular carcinoma and luminal A-like breast cancer, and hence may be driven by sustained hormonal signalling. Molecular analyses reveal a spectrum of factors that could contribute to de novo or acquired resistance to therapy and disease progression.

SDH-deficient renal cell carcinoma associated with biallelic mutation in succinate dehydrogenase A: comprehensive genetic profiling and its relation to therapy response.

Succinate dehydrogenase (SDH)-deficient renal cell carcinoma (RCC) is a rare RCC subtype that is caused by biallelic mutation of one of the four subunits of the SDH complex (SDHA, B, C, and D) and results in inactivation of the SDH enzyme. Here we describe a case of genetically characterized SDH-deficient RCC caused by biallelic (germline plus somatic) SDHA mutations. SDHA pathogenic variants were detected using comprehensive genomic profiling and SDH absence was subsequently confirmed by immunohistochemistry. Very little is known regarding the genomic context of SDH-deficient RCC. Interestingly we found genomic amplifications commonly observed in RCC but there was an absence of additional variants in common cancer driver genes. Prior to genetic testing a PD-1 inhibitor treatment was administered. However, following the genetic results a succession of tyrosine kinase inhibitors were administered as targeted treatment options and we highlight how the genetic results provide a rationale for their effectiveness. We also describe how the genetic results benefited the patient by empowering him to adopt dietary and lifestyle changes in accordance with knowledge of the mechanisms of SDH-related tumorigenesis.

Decreased vancomycin susceptibility in Staphylococcus aureus caused by IS256 tempering of WalKR expression.

Vancomycin-intermediate Staphylococcus aureus (VISA) strains often arise by mutations in the essential two-component regulator walKR; however their impact on walKR function has not been definitively established. Here, we investigated 10 MRSA strains recovered serially after exposure of vancomycin-susceptible S. aureus (VSSA) JKD6009 to simulated human vancomycin dosing regimens (500 mg to 4,000 mg every 12 h) using a 10-day hollow fiber infection model. After continued exposure to the vancomycin regimens, two isolates displayed reduced susceptibility to both vancomycin and daptomycin, developing independent IS256 insertions in the walKR 5' untranslated region (5' UTR). Quantitative reverse transcription-PCR (RT-PCR) revealed a 50% reduction in walKR gene expression in the IS256 mutants compared to the VSSA parent. Green fluorescent protein (GFP) reporter analysis, promoter mapping, and site-directed mutagenesis confirmed these findings and showed that the IS256 insertions had replaced two SigA-like walKR promoters with weaker, hybrid promoters. Removal of IS256 reverted the phenotype to VSSA, showing that reduced expression of WalKR did induce the VISA phenotype. Analysis of selected WalKR-regulated autolysins revealed upregulation of ssaA but no change in expression of sak and sceD in both IS256 mutants. Whole-genome sequencing of the two mutants revealed an additional IS256 insertion within agrC for one mutant, and we confirmed that this mutation abolished agr function. These data provide the first substantial analysis of walKR promoter function and show that prolonged vancomycin exposure can result in VISA through an IS256-mediated reduction in walKR expression; however, the mechanisms by which this occurs remain to be determined.

Comparative analysis of Mycobacterium tuberculosis pe and ppe genes reveals high sequence variation and an apparent absence of selective constraints.

Mycobacterium tuberculosis complex (MTBC) genomes contain 2 large gene families termed pe and ppe. The function of pe/ppe proteins remains enigmatic but studies suggest that they are secreted or cell surface associated and are involved in bacterial virulence. Previous studies have also shown that some pe/ppe genes are polymorphic, a finding that suggests involvement in antigenic variation. Using comparative sequence analysis of 18 publicly available MTBC whole genome sequences, we have performed alignments of 33 pe (excluding pe_pgrs) and 66 ppe genes in order to detect the frequency and nature of genetic variation. This work has been supplemented by whole gene sequencing of 14 pe/ppe (including 5 pe_pgrs) genes in a cohort of 40 diverse and well defined clinical isolates covering all the main lineages of the M. tuberculosis phylogenetic tree. We show that nsSNP's in pe (excluding pgrs) and ppe genes are 3.0 and 3.3 times higher than in non-pe/ppe genes respectively and that numerous other mutation types are also present at a high frequency. It has previously been shown that non-pe/ppe M. tuberculosis genes display a remarkably low level of purifying selection. Here, we also show that compared to these genes those of the pe/ppe families show a further reduction of selection pressure that suggests neutral evolution. This is inconsistent with the positive selection pressure of "classical" antigenic variation. Finally, by analyzing such a large number of genes we were able to detect large differences in mutation type and frequency between both individual genes and gene sub-families. The high variation rates and absence of selective constraints provides valuable insights into potential pe/ppe function. Since pe/ppe proteins are highly antigenic and have been studied as potential vaccine components these results should also prove informative for aspects of M. tuberculosis vaccine design.

SNP/RD typing of Mycobacterium tuberculosis Beijing strains reveals local and worldwide disseminated clonal complexes.

The Beijing strain is one of the most successful genotypes of Mycobacterium tuberculosis worldwide and appears to be highly homogenous according to existing genotyping methods. To type Beijing strains reliably we developed a robust typing scheme using single nucleotide polymorphisms (SNPs) and regions of difference (RDs) derived from whole-genome sequencing data of eight Beijing strains. SNP/RD typing of 259 M. tuberculosis isolates originating from 45 countries worldwide discriminated 27 clonal complexes within the Beijing genotype family. A total of 16 Beijing clonal complexes contained more than one isolate of known origin, of which two clonal complexes were strongly associated with South African origin. The remaining 14 clonal complexes encompassed isolates from different countries. Even highly resolved clonal complexes comprised isolates from distinct geographical sites. Our results suggest that Beijing strains spread globally on multiple occasions and that the tuberculosis epidemic caused by the Beijing genotype is at least partially driven by modern migration patterns. The SNPs and RDs presented in this study will facilitate future molecular epidemiological and phylogenetic studies on Beijing strains.

Evolution of multidrug resistance during Staphylococcus aureus infection involves mutation of the essential two component regulator WalKR.

Antimicrobial resistance in Staphylococcus aureus is a major public health threat, compounded by emergence of strains with resistance to vancomycin and daptomycin, both last line antimicrobials. Here we have performed high throughput DNA sequencing and comparative genomics for five clinical pairs of vancomycin-susceptible (VSSA) and vancomycin-intermediate ST239 S. aureus (VISA); each pair isolated before and after vancomycin treatment failure. These comparisons revealed a frequent pattern of mutation among the VISA strains within the essential walKR two-component regulatory locus involved in control of cell wall metabolism. We then conducted bi-directional allelic exchange experiments in our clinical VSSA and VISA strains and showed that single nucleotide substitutions within either walK or walR lead to co-resistance to vancomycin and daptomycin, and caused the typical cell wall thickening observed in resistant clinical isolates. Ion Torrent genome sequencing confirmed no additional regulatory mutations had been introduced into either the walR or walK VISA mutants during the allelic exchange process. However, two potential compensatory mutations were detected within putative transport genes for the walK mutant. The minimal genetic changes in either walK or walR also attenuated virulence, reduced biofilm formation, and led to consistent transcriptional changes that suggest an important role for this regulator in control of central metabolism. This study highlights the dramatic impacts of single mutations that arise during persistent S. aureus infections and demonstrates the role played by walKR to increase drug resistance, control metabolism and alter the virulence potential of this pathogen.

Rifampicin reduces susceptibility to ofloxacin in rifampicin-resistant Mycobacterium tuberculosis through efflux.

RATIONALE: Central dogma suggests that rifampicin resistance in Mycobacterium tuberculosis develops solely through rpoB gene mutations. OBJECTIVE: To determine whether rifampicin induces efflux pumps activation in rifampicin resistant M. tuberculosis strains thereby defining rifampicin resistance levels and reducing ofloxacin susceptibility. METHODS: Rifampicin and/or ofloxacin minimum inhibitory concentrations (MICs) were determined in rifampicin resistant strains by culture in BACTEC 12B medium. Verapamil and reserpine were included to determine their effect on rifampicin and ofloxacin susceptibility. RT-qPCR was applied to assess expression of efflux pump/transporter genes after rifampicin exposure. To determine whether verapamil could restore susceptibility to first-line drugs, BALB/c mice were infected with a MDR-TB strain and treated with first-line drugs with/without verapamil. MEASUREMENTS AND MAIN FINDINGS: Rifampicin MICs varied independently of rpoB mutation and genetic background. Addition reserpine and verapamil significantly restored rifampicin susceptibility (p = 0.0000). RT-qPCR demonstrated that rifampicin induced differential expression of efflux/transporter genes in MDR-TB isolates. Incubation of rifampicin mono-resistant strains in rifampicin (2 µg/ml) for 7 days induced ofloxacin resistance (MIC > 2 µg/ml) in strains with an rpoB531 mutation. Ofloxacin susceptibility was restored by exposure to efflux pump inhibitors. Studies in BALB/c mice showed that verapamil in combination with first-line drugs significantly reduced pulmonary CFUs after 1 and 2 months treatment (p < 0.05). CONCLUSION: Exposure of rifampicin resistant M. tuberculosis strains to rifampicin can potentially compromise the efficacy of the second-line treatment regimens containing ofloxacin, thereby emphasising the need for rapid diagnostics to guide treatment. Efflux pump inhibitors have the potential to improve the efficacy of anti-tuberculosis drug treatment.

Proteogenomic analysis of polymorphisms and gene annotation divergences in prokaryotes using a clustered mass spectrometry-friendly database.

Precise annotation of genes or open reading frames is still a difficult task that results in divergence even for data generated from the same genomic sequence. This has an impact in further proteomic studies, and also compromises the characterization of clinical isolates with many specific genetic variations that may not be represented in the selected database. We recently developed software called multistrain mass spectrometry prokaryotic database builder (MSMSpdbb) that can merge protein databases from several sources and be applied on any prokaryotic organism, in a proteomic-friendly approach. We generated a database for the Mycobacterium tuberculosis complex (using three strains of Mycobacterium bovis and five of M. tuberculosis), and analyzed data collected from two laboratory strains and two clinical isolates of M. tuberculosis. We identified 2561 proteins, of which 24 were present in M. tuberculosis H37Rv samples, but not annotated in the M. tuberculosis H37Rv genome. We were also able to identify 280 nonsynonymous single amino acid polymorphisms and confirm 367 translational start sites. As a proof of concept we applied the database to whole-genome DNA sequencing data of one of the clinical isolates, which allowed the validation of 116 predicted single amino acid polymorphisms and the annotation of 131 N-terminal start sites. Moreover we identified regions not present in the original M. tuberculosis H37Rv sequence, indicating strain divergence or errors in the reference sequence. In conclusion, we demonstrated the potential of using a merged database to better characterize laboratory or clinical bacterial strains.

Using a label-free proteomics method to identify differentially abundant proteins in closely related hypo- and hypervirulent clinical Mycobacterium tuberculosis Beijing isolates.

Although the genome of the Mycobacterium tuberculosis H37Rv laboratory strain has been available for over 10 years, it is only recently that genomic information from clinical isolates has been used to generate the hypothesis of virulence differences between different strains. In addition, the relationship between strains displaying differing virulence in an epidemiological setting and their behavior in animal models has received little attention. The potential causes for variation in virulence between strains, as determined by differential protein expression, have similarly been a neglected area of investigation. In this study, we used a label-free quantitative proteomics approach to estimate differences in protein abundance between two closely related Beijing genotypes that have been shown to be hyper- and hypovirulent on the basis of both epidemiological and mouse model studies. We were able to identify a total of 1668 proteins from both samples, and protein abundance calculations revealed that 48 proteins were over-represented in the hypovirulent isolate, whereas 53 were over-represented in the hypervirulent. Functional classification of these results shows that molecules of cell wall organization and DNA transcription regulatory proteins may have a critical influence in defining the level of virulence. The reduction in the presence of ESAT-6, other Esx-like proteins, and FbpD (MPT51) in the hypervirulent strain indicates that changes in the repertoire of highly immunogenic proteins can be a defensive process undertaken by the virulent cell. In addition, most of the previously well characterized gene targets related to virulence were found to be similarly expressed in our model. Our data support the use of proteomics as a complementary tool for genomic comparisons to understand the biology of M. tuberculosis virulence.

Multiple, independent, identical IS6110 insertions in Mycobacterium tuberculosis PPE genes.

IS6110 is a transposable element found in Mycobacterium tuberculosis complex members. Regions of preferential IS6110 integration occur within the M. tuberculosis genome but the element has not previously been shown to exhibit any sequence-specific integration preferences. Here we provide evidence for multiple independent IS6110 insertions into identical, or near-identical, positions within the highly homologous 5' region of 3 PPE genes.

Evidence for a rapid rate of molecular evolution at the hypervariable and immunogenic Mycobacterium tuberculosis PPE38 gene region.

BACKGROUND: PPE38 (Rv2352c) is a member of the large PPE gene family of Mycobacterium tuberculosis and related mycobacteria. The function of PPE proteins is unknown but evidence suggests that many are cell-surface associated and recognised by the host immune system. Previous studies targeting other PPE gene members suggest that some display high levels of polymorphism and it is thought that this might represent a means of providing antigenic variation. We have analysed the genetic variability of the PPE38 genomic region on a cohort of M. tuberculosis clinical isolates representing all of the major phylogenetic lineages, along with the ancestral M. tuberculosis complex (MTBC) member M. canettii, and supplemented this with analysis of publicly available whole genome sequences representing additional M. tuberculosis clinical isolates, other MTBC members and non tuberculous mycobacteria (NTM). Where possible we have extended this analysis to include the adjacent plcABC and PPE39/40 genomic regions. RESULTS: We show that the ancestral MTBC PPE38 region comprises 2 homologous PPE genes (PPE38 and PPE71), separated by 2 esat-6 (esx)-like genes and that this structure derives from an esx/esx/PPE duplication in the common ancestor of M. tuberculosis and M. marinum. We also demonstrate that this region of the genome is hypervariable due to frequent IS6110 integration, IS6110-associated recombination, and homologous recombination and gene conversion events between PPE38 and PPE71. These mutations result in combinations of gene deletion, gene truncation and gene disruption in the majority of clinical isolates. These mutations were generally found to be IS6110 strain lineage-specific, although examples of additional within-lineage and even within-cluster mutations were observed. Furthermore, we provide evidence that the published M. tuberculosis H37Rv whole genome sequence is inaccurate regarding this region. CONCLUSION: Our results show that this antigen-encoding region of the M. tuberculosis genome is hypervariable. The observation that numerous different mutations have become fixed within specific lineages demonstrates that this genomic region is undergoing rapid molecular evolution and that further lineage-specific evolutionary expansion and diversification has occurred subsequent to the lineage-defining mutational events. We predict that functional loss of these genes could aid immune evasion. Finally, we also show that the PPE38 region of the published M. tuberculosis H37Rv whole genome sequence is not representative of the ATCC H37Rv reference strain.